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On a Xew Method of Determining the General Perturbations of the Minor 

Planets . 

By W. F. McK\ Ritter , of Milton , Pa. 

(Read before the American Philosophical Society , March 17, 1893.) 

In finding the general perturbations of the minor planets the special 
difficulty arises from the large eccentricity and inclination of these bodies. 
The methods used in case of the major planets fail when applied to the 
minor, on account of want of convergence in the series. Astronomers 
were content, therefore, for a long time, with computing the special per¬ 
turbations of these bodies from epoch to epoch. Hansen finally succeeded 
in effecting a solution of the problem, and his work entitled, Ausein- 
andtrsetzung einer Zweckmassigen Metliode zur Berechnung dev Abso- 
luten Stoungen der Kleiner Planeten, contains all the formula) necessary 
in the cases thus far occurring. 

Instead of determining the perturbations of the coordinates, rectangu¬ 
lar or polar, or of finding the variations of the elements, as had been 
done by his predecessors, Hansen, in his mode of treatment, regards the 
elements as constant, and finds what we may term the perturbation of the 
time. Thus, in place of the time, he uses a function of the time, which 
he designates by z ; so that if g 0 is the mean anomaly at the epoch, we 
have the mean anomaly at any time, in the disturbed orbit, given by 
g 0 -\-n 0 z, ii 0 being the mean daily motion, and being one of the constants. 
If there were no perturbations we should have g Q -\-n^t, t being the time 
elapsed since the epoch. 

In effecting his solution of the problem, Hansen does not attempt to 
give general and complete analytical expressions of the series. Instead, 
he, at the start, converts the coefficients into numbers, and multiplies the 
series together, two and two, by the methods of trigonometry. Thus, 
although we find, finally, the perturbations as functions of the time, that 
is, have the general perturbations, yet, in applying the method to differ¬ 
ent bodies, we must find the values of all the quantities involved for the 
particular case under consideration. It would be a great advance if wc 
had at hand complete analytical expressions, of sufficient convergence, as 
is the case with the larger planets. 

Besides the method of multiplying series together by the methods of 
trigonometry, which Hansen calls “Mechanical Multiplication*'—a 
method he was the first to employ—he also adopts different angles with 
which to express his arguments. Thus at the outstart he uses the eccen¬ 
tric anomaly for both bodies. When he has computed the powers of the 
reciprocal of the distance between the disturbed and disturbing bodies, 
lie transforms from eccentric to mean anomaly in ense of the disturbing 
body. And then, when lie lias expressions for the perturbing function 
and the forces, he makes another transformation so as to be able to effect 
the Integrations. 
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The transformations must be clone with great care, and require a large 
measure of time. In addition to the tedium arising from extended opera¬ 
tions of this kind (which must be generally done in duplicate to insure 
accuracy), many of the processes in various stages of the w T ork are not 
easily grasped, and certainty is often only secured by performing the 
numeiical calculations. Thus, then, although the method has been pub¬ 
lished for a long time, it has been applied only in a very limited number 
of cases. Watson, in the Preface to his Theoretical Astronomy, says : 
“The refined and difficult analysis and the laborious calculations involved 
were such that, even after Hansen’s methods were made known, astron¬ 
omers still adhered to the method of special perturbations by the variation 
of constants as developed by Lagrange.’* 

Hansen seems himself to have felt the force of these drawbacks on his 
method, as in a posthumous memoir devoted to the larger planets he 
abandons his peculiar method of treatment and uses that of Lagrange. 

As far as the minor planets are concerned, there is no doubt that Hansen’s 
method, as left by him, is too long and difficult to be practicable. 

What we need now is some mode of determining general perturbations 
that is easily applied and sufficiently short to attract the efforts of a larger 
number of competent computers. Only in this way can the constantly 
growing material be utilized. The new method of treatment will now be 
given as briefly as possible. 

If A be the distance betwceu the disturbed and disturbing bodies, 
Hansen has the equation 

n _ 1 

(j) ~ { °~ q cos — } 2 { i— ?i( cose '+ £>} 2 

forfinding (*), (~) , etc. 

Instead of the two factors of the second member, I have used a trans¬ 
formation of them given by Hill, and have 
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The coefficients of both these factors are the La Place coefficients, and 
their values have been tabulated. Thus the part of the work relating to 

the determination of expressions for c> O' , etc., is rendered 


comparatively short and simple. 

In finding A in terms of the radii vectores of the two bodies and of the 
cosine of the angle between these radii-vectores, the true anomaly of both 
bodies is introduced. In the analysis we use the equivalent functions of 
the eccentric anomaly for those of the true anomaly, and then, when 
making the numerical computations, we cause the eccentric anomaly of 
the disturbed body to disappear. This is accomplished by dividing the 
circumference into a certain number of equal parts relative to the mean 
anomaly, and employing for the eccentric anomaly its numerical values 
corresponding to the various values of the mean anomaly. 

Having found the expressions of f“Y (1Y , etc., in series, in which 


the angles are the mean anomaly of the disturbed and the eccentric 
anomaly of the disturbing body, the series are changed at once into others 
in which both angles are mean anomalies. To effect this transformation 
there is need of functions called the J functions ; and a chapter is given in 
which the expressions for these functions arc found in a lorm convenient 
for application. 

When we have the powers of the reciprocal of the distance between 
the disturbed and disturbing bodies, we next find the term expressing the 
effect of the action of the disturbing body on the sun. This is effected 
without difficulty. 

The expressions for the perturbing function and the perturbing forces 
can now be formed. Instead of using the force involving the true anomaly, 
the transformation of this, in which the mean anomaly appears instead of 
the true, has been used. This is the method givcu by Hansen in his post¬ 
humous memoir, in which he has abandoned some of his former notions. 
The disturbing forces employed arc those in the direction of the disturbed 
radius-vector, in the direction perpendicular to this radius-vector, and in 
the direction perpendicular to the plane of the orbit. The forces in these 
three directions have been deduced from those in the direction of the 

three rectangular axes. The force a is found at once from the per- 

dg 

turbing function by difierentiating with respect to the mean anomaly, g. 


To find the other two forces symbolized by a r. and « 9 z being 

dr dz 

the coordinate perpendicular to the plane of the orbit, it is necessary to 
multiply a number of series together, two and two, by the formula! of 
plane trigonometry. 

Having the values of the forces, we next find the value of a function 
W obtained by the integration of the expression 

• lW ■ J. n n , d & 
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A and B being two factors easily determined. W being known, the func¬ 
tion W is next found by simple mechanical processes, and the perturba¬ 
tions of z and of the radius-vector are found at once by the equations 


dt 

rdW 


n. tiz = n J~. W. 

v = -l u J-dr- dt ’ 

Y being a particular form for g . 

The symbol d designates the perturbation of the quantity to which it is 
prefixed. 

The perturbation of the latitude is found by integrating the equation 

d.JL 

_ cos i- n a , d 3, 

n. dt ” a dz 

where C is a factor found in the same manner that A and B were. 

To find n.8z, or the variation of the mean anomaly, two integrations 
are necessary ; in finding the perturbation of the radius-vector, and of the 
latitude, one integration is needed for each. 

The arbitrary constants introduced by these integrations are so deter¬ 
mined that the perturbations become zero for the epoch of the elements. 

In making an application of his formulae, Hansen selected the planet 
Egeria, whose eccentricity is comparatively small, the angle of eccen¬ 
tricity being less than five degrees. In making use of the formulae given 
in the method here presented the eccentricity is considerably larger. The 
convergence of the series is, however, all that can be desired. In com¬ 
puting the perturbations of those of the minor planets whose eccentrici¬ 
ties and inclinations are quite large, it may be necessary to divide the 
circumference into a larger number of parts. In exceptional cases, such 
as for Pallas, it may be necessary to divide the circumference into thirty-two 
parts. In case of the applications made of the present method, sixteen 
divisions have been used ; this is the number employed by Hansen for 
Egeria. 

When a larger number than sixteen is used, the calculation of the 
values of O (1)' , etc., is longer; the process is, however, the 


same in every case. 

After the perturbations have been found it is necessary to have them 
in convenient form for the computation of epbemerides, and there has, 
hence, been added the method employed for doing this. 

The writer has endeavored to present the whole theory in a manner 
easily comprehended by those having a respectable mathematical educa¬ 
tion, and in a compass such that the computations can be performed 
within a reasonably short time. The endeavor throughout has been to 
use convenient methods, not to devise new ones. 


